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Forced oscillation of a system composed of two pendulums coupled by a spring in the presence 
of damping is investigated. In the steady state and within the small angle approximation we solve 
the system equations of motion and obtain the amplitudes and phases of in terms of the frequency 
of the sinusoidal driving force. The resonance frequencies are obtained and the amplitude ratio 
is discussed in details. Contrary to a single oscillator, in this two-degree of freedom system four 
resonant frequencies, which are close to mode frequencies, appear. Within the pass-band interval 
the system is shown to exhibit a rich and complicated behaviour. It is shown that damping crucially 
affects the system properties. Under certain circumstances, the amplitude of the oscillator which is 
directly connected to the driving force becomes smaller than the one far from it. Particularly we 
show the existence of a driving frequency at which the connected oscillator's amplitude goes zero. 



I. INTRODUCTION 

Forced oscillation of many body systems constitutes 
an important and interesting subject in wave phenom- 
ena. This subject is normally discussed incompletely in 
the text books and various aspects have remained unex- 
plored mill. In reference [J this issue is discussed 
in the simplest framework of a system having two de- 
grees of freedom consisting of two pendulums connected 
to each other by a spring. This system has two mode fre- 
quencies uji and UJ2 > lui. It is qualitatively argued that 
if the driving frequency is lower that uji or greater than 
UJ2 the amplitude of the oscillator directly connected to 
the driving force is larger than the one far from it. The 
further the driving frequency from the mode frequency, 
the larger the amplitude ratio. In this manner the system 
behaves as a mechanical filter which allows frequencies in 
the pass band region [uji , UJ2] ■ Reference [1] does not dis- 
cuss what occurs in the interesting case where the driving 
frequency lies in pass-band region and only states that in 
the pass-band region the amplitudes are comparable to 
each other. In this paper we aim to quantitatively inves- 
tigate the properties of such system in the entire range 
of e [0, oo[. We show interesting features will arise 
when Lo lies in the pass-band region. In particular in the 
absence of damping it is shown there exists an interme- 
diate frequency at which the amplitude of the connected 
oscillator to the driving force becomes zero. 



II. 



STATEMENT OF THE PROBLEM 



Consider two pendulums a and b each having length 
L and bob mass M which are coupled via a spring with 
constant K as illustrated in figure (1). Pendulums are 
subjected to frictional damping forces and each one has 
the same damping constant per mass F. Damping can 
be provided by something rubbing on supporting strings 
or the bobs. Imagine we derive the system by attaching 
one of the pendulums say a to a sinusoidal driving force 
with amplitude per unit length Fq and frequency w. The 
equations of motion can easily be written in small angle 




Pendulum a 



Pendulum b 



FIG. 1: Two coupled pendulums with masses M which are 
derived by one end by a sinusoidal force. 



approximation: 



A/*, = -^*,-X(^,-*b)-Mr^, + FoCoswt (1) 



Mg 



(2) 



In which ^'q and are the angular displacements of 
pendulums a and b. Introducing the normal coordinates 
+ and *2 = 5(*a - *a) and normal 



^1^1 — 

*i 2 

frequencies uJi — y o,iiu ^/^^ — , / -j- i -jj- 
equations (1) and (2) into the following decoupled form: 



y- and ll)2 



2K 



we recast 



2M 



cos Ojt 



^2 + + UJo'^2 = cos 

2-r 2 2 2M 



(3) 



(4) 



In steady state when the effect of transients are negli- 
gible and ^'2 appear to be as two independent forced 
harmonic oscillators with frequency w: 
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^'i = Ai{uj) cos{iot - Si{lu)) 

^'2 — A2{uj) cos{ujt — (^2(0;)) 
Amplitudes are obtained as follows [Ij: 



A,{lo) 



Fo 



(5) 
(6) 

*=1,2 (7) 



Trigonometric functions of the phase differences re- 
spect to driving force turn out to be: 



sin Si (cu) 



cos 5i{uj) 



tan 6i{uj) 



v/K'-Cj2)2+r2w2 



* = 1,2 (9) 



i = 1,2 



(10) 



The steady state displacements of oscillators a and b 
will be = '^1 + ^2 and = *i - ^1/2 respectively. 
After some straightforward trigonometric manipulations 
we obtain: 



*a = Aa COa{ujt - Sa) 

In which the amplitude Aa and phase Sa are: 



(11) 



Aa^ jAl+Al + 2AiA2Cos{Si-S2) (12) 



Ai sin (5i + ^2 sin S2 

smda = 
cos (5a = 
tan (5a = 

In a similar fashion we obtain for oscillator b: 

^'b = Aij cos{ujt — Sb) 





Aa 




Ai 


cos(5i + A2 


cos (52 




Aa 




Ai 


sin(5i + yl2 


sin S2 


Ax 


cos S\ + yl2 


cos (52 



(13) 



(14) 



(15) 



(16) 



sin Sb 



cosSb — 



tan Sb 



Ai sin Si — A2 sin S2 





Ah 




Ai 


cos(5i — A2 


cos S2 




Ab 




Ai 


sin(5i — A2 


sin (52 


Ai 


cos(5i — A2 


cos (52 



(18) 



(19) 



(20) 



The main questions we would like to address are as 
follows: 



i = l,2 (8) frequency a; 



1) What is the dependence of ^ 



on the driving 



2) What are the dependence of phases Sa and Sb on 
the driving frequency uj ? 

We address these questions in two separate cases. In 
the first case we assume the drag force is absent i.e.; 
r = 0. In the second case which is far more challenging 
the drag force is non zero T ^ 0. Let us first discuss case 
one. 



III. NO DAMPING FORCE 

In this case where F = the mode amplitudes and 
phases are simplified as follows: 



A,. 



Fo 



2M|a;2_^2| 



; sin Si — tan (5; = (21) 



cos (5,- — 



|W? - C^2 



1,2 



(22) 



Three regimes are normally discussed in the literature: 
Lu < uji < ijj2] uJi < UJ < UJ2 and uji < uj2 < lj. We discuss 
each regime separately. 



A. Regime I: < cji < tJ2 

In this regime where w is below the pass-band lower 
cut-off we have: cos(5i = cos (52 = 1 which give Si = S2 = 
0. By (13) we conclude sin (5a = which implies Sa is 
either zero or tt. The correct choice depends on cos (5a. 
We have cos Sa = ^^j^- On the other hand Si = S2 = 
and (12) give Aa = Ai + A2 which in turn gives rise 
to cos (5a = 1. This excludes Sa = tt and therefore we 
conclude Sa = 0- Similarly we have sin Sb = which gives 
Sb is either zero or tt. By (19) we have cos Sb — 



Ab= JAI+Al- 2A1A2 cos((5i - S2) (17) On the other hand (17) gives A = \Ai - A2\. Having in 
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mind that in regime I we have lo < uJi < UJ2 together with 
(21) we can conclude that Af, = Ai — A2 which gives rise 
to cos^b = 1. This fixes Sb = 0. We therefore conclude 
that in regime I both oscillators a and b are in phase with 
the driving force. Lastly we evaluate the amplitude ratio 
According to our findings we have 



A, 
At 



A1+A2 
A1-A2 



> 1 



(23) 



This means oscillator a which is directly connected to 
driving force has a larger amplitude. 



B. Regime II: loi <u) < L02 

In this regime we can simply determine 6b and Ab- 
According to (22) cos^i = — 1 and cos(52 = 1 which lead 
to Ab — Ai + A2. By (19) we have cos (5b = —1 which 
gives Sb = TT. However, determination of 6 a and Aa is 
not so simple as that was for Sb- We note cos^i = — 1 
and cos52 = 1 give Si = w and <52 = 0. Moreover, 
sinJo = which gives Sa — or n and cosSa = ^E^- 
To determine the sign of the denominator we note 

that according to (21) if w < \f^^ then Ai > A2 
which fixes cos Jo = — 1. This yields Sa = tt. Since we 

will frequently encounter the frequency y ^ 
remainder of the paper from now on we denote it by w. 
On the other hand if w > w we have Ai < A2 which fixes 
cos^a = 1. This yields Sa = 0. In the case of equality 
i.e.; LJ = ui we have Ai = A2. 

To summarize, when the driving frequency lies in 
pass band that is between the mode frequencies three 
distinctive cases are identified: 

Case 1) ui < u < ui 

In this case Sa = tt (oscillator connected to the driving 
force is out of phase with it !), Sb = n, Aa = \Ai — = 
Ai — A2, Ab = Ai + A2. Therefore the amplitude ratio 
reads: 



Aa ^ A1-A2 

Ab Ai + A2 
Case 2) ui <oj < 0J2 



< 1. 



(24) 



In this case Sa = and Aa = \Ai — A2I = A2 — Ai, 
Ab — Ai + A2. Therefore the amplitude ratio turns out 
to be: 



Aa ^ A2-A1 

Ab Ai + A2 
Case 3) ui = oj 



< 1. 



(25) 




FIG. 2: Dependence of ratio ^ on a; in absence of damping 
for three values of coupling parameter K. 



In this case Aa = that is the connected oscillator to 
the driving force remains immobile, its phase Sa is there- 
fore undetermined. This case is interesting and intro- 
duces a frequency w at which the amplitude of the oscil- 
lator connected to the driving force is zero while the sec- 
ond oscillator has a non zero amplitude Ab = 2Ai = 
Therefore the amplitude ratio ^ is zero. 



C. Regime III: wi < a;2 < w 

In this regime the driving frequency is above the pass 
band upper cut off. We simply obtain cosJi = cos52 = 
— 1 which give Si = 82 = n. By a similar argument we 
obtain Sa = n, Sb = 0, Aa = Ai + A2 and Ab = A2 — Ai. 
The amplitude ratio becomes: 



Ab 



A1 + A2 
A2-A1 



> 1. 



(26) 



Note that in this case oscillator a is out of phase 
with the driving force. This aspect is discussed in 
section 3.3 of reference [1] which qualitatively predicts 
that when the system is derived above the high cut-off 
frequency UJ2 its constituents phases resemble to that 
of the highest mode of free oscillation (out of phase 
relative to each other). Figure (2) shows the dependence 
of amplitude ratio on u) for various values of coupling K. 
Throughout the paper we have used the unit in which 
M = g = L = 1. In this unit wi turns out to be one. 
Moreover, we have set K = 10 unless otherwise stated. 
For this value of K we have U2 = 4.58. 
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IV. EFFECT OF DAMPING FORCE 



A. Regime I: u < oji < L02 



Now we turn to the more realistic case where damping 
forces are present i.e.; F 7^ 0. By substitution of (8) and 
(9) into (15) and (20) we obtain: 



tan^a = 



tan^ft = 



[(a;2_^2)(^2_^2)_r2a;2]- 



(27) 



(28) 



Four independent polynomials appear in the above ex- 
pressions. Let us name them as follows: 



q{uj) = (wf - co'){cji - cj') + r'uj\ 



p{uj) = (Wi - W^)(W2 - W^) 



■ rV. 



(30) 



According to tables (1-7) when oj is below the lower 
cut-off uji, q and h arc positive so wc conclude that 
tan Sa > {5a > 0) which means the oscillator connected 
to driving force is in phase with it. To determine the sign 
of p, we introduce two frequencies 



.(p) 



■ \— h 

^ 2 ^ 



In which A^^^ is defined as follows: 



A(p) = (w2^a;i+r2)2-4a;2 



,2. .2 



(35) 



(36) 



(p) 

In appendix we show w_ < wi and p > for 

< w < Lo^^ whereas p < for cj^^ < cj < cji . 

Therefore if < w < then tanJ^ > > 0. 

(29) If wi''^ < Lu < LUi then UiiiSb < =^ Sb < 0. We note 

that 00^ > cJi therefore it lies outside the region under 
consideration. Now let us discuss the amplitude ratio 
4^ . As shown in tables (3-7) three regimes of T should 



be distinguished: 



S{w) = (wf - CO^y + (^2 

With these definitions, we have: 



.2\2 



-I- 2r^a; 



2, .2 



- Tujs - Fujh 

tan da = — ;— tandb = . 

qh p 



(31) 



(32) 



(33) 



Let us first discuss amplitudes ratio. According to (12) 
and (17) ^ depends on cos(5i — 62)- If cos((5i — ^2) > 
then Aa > otherwise if cos(5i — ^2) < then Aa < Af,. 
Using (8) and (9) we obtain after somewhat algebra: 



cos((5i - 82) 



(34) 



Since the denominator of (34) is positive, the sign of 
cos((5i — ^2) is the same as q. The signs of 6a and 8^ 
depends on signs of p, g, s and h. Clearly s > so we 
only have to determine the signs of p. q and h in the range 
< w < 00. Polynomial h is quadratic in w therefore it 
would be easy to determine its sign. By contrast, p and 
q arc quartic and their sign determination are not trivial. 
This task is done in details in appendix and the results 
are shown in tables (1-7). We only quote and use the 
results here. Analogous to the case of free damping case 
we discuss three regimes: 



A: < r < a;2 - wi 

Here have > which implies cos(5i — ^2) > 
Therefore the amplitude ratio ^ is greater than one. 

B : r = a;2 — wi 

Here have g < therefore cos(5i — ^2) < and 
therefore the amplitude ratio ^ is less than one. 

C : r > 0)2 — Wi 

Here g > therefore cos((5i — 82) > ^ and hence the 
amplitude ratio ^ is greater than one. 



B. Regime II: oji < w < W2 

According to Appendix, sign determination of q cru- 
cially depends on F. Let us first discuss tan (5^. According 
to (33) the sign of tan 8\, is identical to ^ . For wi < w < w 
we have tan^b < whereas for w < a; < 0^2 tan (5b be- 
comes positive. Now we discuss tan 8a- According to (33) 
the sign of tan 8a is the same as qh. We introduce two 
frequencies 



? + a;2 _ p2 ± 



, ,(9) _ 



I ' which a'') is defined as follows: 



(37) 
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A(9) = (a;?+a;i-r 



2\2 



(38) 



Five states are identified: 
A: < r < -^^L= 

\/2(w?+wi) 

According to appendix we have tan (5a > for 
u>i < u> < u)^^ and w < w < w^j?' and tan (5a < in 
the remainder of the interval [a;i,a;2]- Now we discuss 



ampHtude ratio. In the interval [u;i,a;i''^[ and ]a;^'',a;2] 
we have q > which gives 4^ > 1. In the remainder of 



,(9) 



the interval ^ < 1. At w = Wj?^ we have ^ = 1 

2 2 

B : r = 



\/2K+wi) 

According to appendix we have tan^^ > for 

LOi < Lo < wL'"* and tan (5a < in the remainder of the in- 
terval [cji; CJ2]- The amplitude ratio is identical to part A. 



C : 



In this range of T tan (5a > for wi < w < u>^^ and 



,(<?) 



< uj < Co and tan (5a < in the remainder of the 



interval [ci;i,u;2]- The amplitude ratio is identical to that 
in part A. 



D : r 



: W2 — Wl 



At this particular value of F tan (5a > for two sub 
intervals: Wi < w < w* and uj* < oj < u). tan (5a < in 
the remainder of the interval [a;i,a;2]. Note that in this 
case, at two special frequencies u* = y^tJIZJi and u) 
the expression qh becomes zero and hence tan (5a goes 
infinity. The notable point is the amplitude ratio ^ 
is greater than one except at ui* = y/uiiUJ2 where it 
becomes unity. 

E : T > ui2 — uJi 

In this range of F tan Sa > for wi < cj < cj and 
tan (5a < in the remainder of the interval [wi, ^2]- The 
amplitude ratio ^ is greater than one in the whole range 

[uil,UJ2]- 



C. Regime III: wi < W2 < w 

We now discuss regime III which is much simpler than 
regime II. In regime III the driving frequency is above the 
upper cut-off frequency W2 . In this regime q > and h < 
0. Therefore (33) gives Sa < 0. According to table (2), 



tan (5;, > for ix)2 < cj < uj^_^^ and tan St < for w'^^ < uj. 
Note that at oj = w'j^'' we have tan 6h = ±00 depending on 
how w approaches w^"* . The amplitude ratio ^ is always 



(p) 




FIG. 3: Dependence of ratio ^ on in presence of damping 
for some values of damping constant F in states A to E. 



greater than unity for W2 < w according to appendix 
tables. 



D. Results and discussion 

In this section we discuss about the role of damping 
and drive on the oscillators characteristics. By substitu- 
tion of the modes amplitudes and phases into (12) and 
(17) we arrive at the following relations for the oscillator 
amplitudes Aa and Ai,: 



Aa 




(w2 _ ^2)2 + r2w2)^^2 _ ^2)2 + ^2^2) 



(39) 



Fo / + w| - 2wf 

2M Y (wf _ £^2)2 + r2a;2)(a;2 - a;2)2 + r2a;2) 



(40) 

To find the amplitude resonant frequencies we should 
set the derivatives of Aa and Ab respect to ui equal to 
zero. For Aa this leads to a polynomial equation of order 
10. For Ai, this task leads to a cubic equation for uj with 
lengthy coefficients. Despite availability of analytical 
solution for a cubic equation, we prefer to proceed 
numerically. Figure (3) depicts the amplitude ratio in 
the whole range of oj for various values of F associated 
to the states A, B, C , D and E described in regime II. 
It is worthwhile mentioning that at F = uj2 — ijJi the 
amplitude ratio is an increasing function of lo whereas in 
other values of F there is minimum value for the ratio. 

Figure (4) shows the dependence of amplitudes Aa and 
Ab on w for a small damping constant. There exist two 

peaks at which both amplitudes arc maximised. The 
peak values of w are resonant frequencies and are very 
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FIG. 4: Dependence of amplitudes Aa and Ai, on a; for a 
small value of damping constant F = 0.3. 



FIG. 6: Dependence of taxiSb on w for two values of damping 
constant F. 




close to (jji and uj2. By increasing F the resonant fre- 
quencies of oscillators a and b begin to deviate from each 
other. In fact, each oscillator posses its own resonant 
frequencies. To sec this more explicitly in figure (5) we 
have drawn amplitude dependence on w for larger values 
of r. Separation of resonant frequencies is apparent. 

Figure (6) exhibits tan^b in the whole range of uj for 
two values of F. 

Notice that tan^b exhibits two jumps (discontinuities) 
for each value of F. These jumps occur at roots of p{ijo) 

which are u:^ . In figures (7-11) we have shown taxi 5a vs 
(jj for those values of F in figure (3). 

In contrast to t&nS), the number of jumps in tan^a 
depends on the value of F. We can have one, two and 
three jumps in figures (7-11). We always have a jump 
at w = (D which is the root of h{u)). Other jumps are 
associated to w^?'' (roots of g(w)) which their existence 
depend on the value of F. 



tgSa 
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FIG. 8: Dependence of tanJa on w for damping constant 
F = 3.015 in state B. 
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VI. APPENDIX 

In this appendix we determine the sign of polynomials 

h,p and q in variation of w. Since h is quadratic in lo 
its sign can simply be determined. For later purposes 



we need to specify the location of u) 



respect 



to Wi and W2- According to the following inequality we 
arrive at the following sign table for h: 



FIG. 9: Dependence of tan 5a on w for damping constant 
r = 3.2 in state C. 



(41) 



tg5a 
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FIG. 10: Dependence of tanJa on w for damping constant 
r = 3.58 in state D. 



+ 



+ 



0J2 



(42) 



A. Sign determination of p 
Introducing f2 = polynomial p can be written as: 



(43) 



Eq. (43) is now quadratic in il. Its discriminant A^^*) 
turns out to be: A^p) = (a;^+w|+r2)2-4w?w|. Utilizing 
the inequality (wi — W2)^ + F^ > wc can simply show 
that A^*') > 0. This implies that equation (43) posses 
two real solutions: 
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FIG. 11: Dependence of tani^o on uj for damping constant 
F = 5.5 in state E. 



^ 2 



(44) 



On the other hand we have : ujI +^2+^^ > V A^p^i 



which gives rise to 
two frequencies: 



> 0. Therefore we can introduce 



(p) _ 



•,(p) 



(45) 



Now the sign of p is easy to determine: p is negative 



if LJ. 



< CO < OJ. 



and positive elsewhere. At uj^^ poly- 



nomial p becomes zero. Next wc should determine the 
positions of loi and W2 respect to uj^^ . We claim: 



W_ < (jJi. 



(46) 
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The proof is quite simple: squaring both sides of (46) 
and using the definition of cj'f' ^ we arrive at the inequal- 
ity (w| — ujf + r^) < \/A- Squaring both sides of this 
inequality after some straightforward algebra we reach 
to a true statement > which completes the proof. 
Similarly we can show that: 



r < ^2 - wi (53) 

If + w| - < i.e., if ^/uJ'f + ijjI < T condition 
(52) reduces to 



iof>uj^. (47) 
In a rather similar fashion we can prove: 

w^^^ > (48) 

The following table depicts the sign determination of 
polynomial p. For the sign of 5^ we have to determine 
the sign of | hence we simultaneously exhibit the signs 

of h and - in the following table. 



T >0J2+uJi (54) 

B) A(«) = 0: 

In this case we have a degenerate real root: 

0- = ^^i^ (55) 

Putting A^'^) equal to zero we should solve the following 
equation to find out the condition. 



h 


+ 


+ 


+ 








P 


+ 










+ 


h/p 


+ 






+ 


+ 





(49) 



B. Sign determination of q 

Introducing = cj^ polynomial q{u)) can be cast into 
the following form: 

q = if - {oj^+u^- r2)f2 + cj^wi (50) 

In terms of ft polynomial q becomes quadratic with 
discriminant A^''^ = {cjf + — F^)^ — AloIljI- Contrary 
to polynomial p where its discriminant was shown to be 
positive, here A^'^ can take any sign which complicates 
the problem. Therefore wc should determine the sign of 
in priori. Three distinctive cases arise: 

A) A(9) > 0: 

In this case we have two real roots: 

Let us find out the condition under which A^^^ > 0. 
Using the definition of A^^^ it can be easily shown that 
the suHicient condition for positivity of A^^^ is: 

+w| - F^l > 2wiW2 (52) 

If + - F^ > i.e.; if < F < ^JiJl + lj^ condition 
(52) reduces to 



|wf + -F^l = 2wiW2 (56) 

If + w| - F^ > (0 < F < ^/ij'f+uj'i) (56) gives 
F = W2-W1 whereas if Wi-|-w|-F^ < {^/iof+^ < F) 
we have F = W2 + wi- 

C) A(«) < 0: 

In this case (50) has no real root. When A^^^ < we 

have \ojf + cj| — F^| < 2cjicj2- It can simply be shown 
that this condition is equivalent to W2— wi < F < lo2+'-^i- 

Now we can determine the sign of A^^) in the entire 
range of < F < 00. The following table exhibits the 

sign of A'-'^h In this table we also show the sign of the 
expression /(F) = ojf +CO2 —T^ which appears in (51). 



F 


UJ2 — COl 


V2u} 


OJi + OJ2 


A(9) 


+ 




+ 


/(r) 


+ + 







(57) 



Note that in the interval F €]a;2 — wi , a;2 + wi [ within 
which A*^') < the sign of q{uj) coincides with the sign of 
the highest coeHicient in (50) (positive). According to the 
above table for < F < W2 — wi we have A^') > hence 
both fl'f^ and 0^5^ are positive (by the very definition of 
A^*') and we can define: 

= ^/^. (58) 

By contrast, in the interval F > 0^2 + both Sl^' 

and fi!^' are negative. This means that forth order 
polynomial q (in terms of oj) has no real root therefore it 
is either positive or negative. Since the coefficient of the 
highest term is positive we conclude that q is positive in 
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the entire range < w < oo. The points T = uj2 — uii 
and T = u)2 + uii need special attention. At F = a;2 — Wi 

we have fl* = — — ^ = u!-i_u!2 > therefore we 

conclude that q{uj) = {uj'^ - fi*)^ = (w^ — uJiUJ2)'^ > 0. 
At the second special point F = a;2 + wi we 



have n* 



q{u)) = (w^ +ujiu)2) > 



= —0J1LO2 < therefore 



The next stage is to determine the locations of uii 
and u)2 respect to u)^^ and lo* = = ^wiW2- It is 
sufficient to consider two cases: < F < W2 — Wi and 
r = a;2 — wi since in the remaining range of F there is 
no real root u for q. 



A: < F < W2 - wi 



In this range of F we claim ui'^^ > uii. To show this we 

square both sides and use the definition of lj'^^ in (45) to 
reach this inequality: 



F^ < 



(59) 



It is not difficult to show that the left hand side of (59) 
is positive for the range < T < LJ2 — t^i- Den oting the 
left hand side by 7 we have 7 > for < F < \/a;f — . 
On the other hand wc have uj2 — ^1 < \/^2 " ^1 there- 
fore 7 > in the range < F < W2 — wi- Now we can 
square both sides of (59) without reversing the inequality 
sign. After some straight forward algebra we arrive to 
a true statement F^ > which proves our claim. In a 
similar fashion we can show that co^^ < C02- Now we turn 



to co^^\ Since 



< u!^^ and > lji we conclude 



that > oJi. We claim that < W2- The ar- 

gument is similar to those for wf?^ and we omit the proof. 

B: r = 0J2 — oJi 

Substituting this value of F into (55) gives fi* = UJ1OJ2 
which yields uj* ~ ^ix)\lo2- This means w* is the geo- 
metrical average of mode frequencies therefore we have 
wi < w* < L02- Subsequent to our previous discussions, 
for F > 0^2 — tt'i we have g' > 0. 

We are not yet at the end of the story. The point is that 
the values of w at which (7, -p and h change sign do not 
coincide. This makes the problem complicated, we recall 
that to find the signs of 8a and bh we have to determine 
the sign of qh and - respectively. For this purpose, we 



should determine the location of w = y ^ ^ ^ respect to 

w2'\wj?^ and w* correspondingly. Now we discuss the 
sign of qh which gives the sign of tg^a- As there are 
three regimes for F we discuss them separately. In the 
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first regime i.e.; < F < W2 — i^i we claim lo 
This is evident by the definition of wi*'' as in (51,58). To 



determine the location of w respect to we rewrite 
u}-^ as follows which is more suitable for comparison: 



^A(9) - F2 



(60) 



Clearly if ^jW) > F^ wc have uj^ > ui and vice versa. 

(9) _ r, 



In the case of equality i.e.; V A(9) 
which gives (by definition of A(«)) 



F^ then 



The condition VAW > F^ impUes: 



(61) 



F< 



(62) 



Similarly the condition V A(«) < F^ implies : 



F > 



(63) 



It is straightforward to show 



< 0^2 — . 



At the specific point F = W2 — ^^i wc can easily show 
that w* = yjLUiui2 < w. We are now able to exhibit the 
simultaneous sign tables of q, h and hence tgSa- 

I: 0<F< 



CO 









(9) 

CO CO^' CO2 


q 


+ 


+ 




+ + 


h 


+ 


+ 


+ 




{qh)-^ 


+ 


+ 




+ 










(64) 


II: F = 


\/2(w?-Fwi) 








CO 





UOi 




(9) 

Ul = CO^ LU2 


q 


+ 


+ 




+ + 


h 


+ 


+ 


+ 




{qh)-^ 


+ 


+ 







(65) 



III: 



CO 







<r < C02 — ooi 



COi 



CO 



(9) 



00. 



(9) 



CO 



L02 



Q 


+ 


+ 




+ 


+ + 


h 


+ 


+ 


+ 


+ 





(66) 
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IV: T = UJ2 — iOi 



U) Ul UJ* UJ IjJ2 

q + + + + + 

h + + + - -_ 

Wi}^ + + + - ~ 

(67) 

V: r > a;2 — Wi 
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+ 


+ 


+ 
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+ 


+ 








+ 
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